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1.      INTRODUCTION 

This  report  is  prepared  in  accordance  with  U.S.  Department 
of  Transportation  contract  DOT-FH-11-7769 .   Under  this  contract, 
Systems  Technology  Associates,  Inc.  provides  field  calibration 
services  for  a  data  acquisition  system  and  reduces  data  from  the 
system  to  investigate  wind  velocities,  directions  and  effects  at 
long-span  bridge  sites.   The  data  acquisition  system  has  been 
installed  on  the  Narragansett  Narrows  bridge  in  Newport,  Rhode 
Island,  and  on  the  Sitka  Harbor  bridge  in  Sitka,  Alaska.   The 
geographic  locations  and  photographs  of  these  sites  are  shown  in 
Figures  1  through  3.   Selected  data  recorded  at  Newport  and 
Sitka  are  presented  in  this  report. 

Section  2  presents  wind  velocity  spectra  from  data  record- 
ed on  four  separate  occasions  between  1971  and  1976.   Data  were 
selected  to  show  the  most  significant  wind  conditions  at  the  two 
bridge  sites  with  respect  to  mean  wind  velocity,  data  variance, 
and  wind  angle.   Section  3  presents  bridge  motion  data  in  the 
form  of  acceleration  spectra.   These  show  predominant  frequen- 
cies of  bridge  motions  during  a  single  event  recorded  at  the 
Sitka  site.   No  bridge  motion  data  was  recorded  at  the  Newport 
site  due  to  instrumentation  difficulties. 

Section  4  contains  a  summary  of  activities  during  the  re- 
porting period  and  describes  future  project  plans.   The  appen- 
dices include  a  chronology  of  significant  project  events  since 
the  inception  of  the  project,  a  summary  of  data  collected  and 
maintained  during  the  project,  and  detailed  data  associated 
with  Sections  2  and  3.   The  remaining  paragraphs  of  this  sec- 
tion briefly  introduce  the  data  acquisition  system  and  the 
computer  programs  used  to  process  project  data. 

1.1    DATA  ACQUISITION  SYSTEM 

Data  is  acquired  using  instrumentation  and  a  special- 
purpose  data  recording  system  which  can  be  installed  in  various 
configurations  at  selected  bridge  sites.   Figures  4  and  5  show 
the  installation  configurations  at  the  Newport  and  Sitka  sites. 

Instrumentation  basically  consists  of  anemometers  to 
measure  wind  velocities  and  directions,  and  accelerometers  to 
measure  bridge  motions.   Anemometers  measure  three  orthogonal 
wind  component  velocities  -  along,  across  and  vertical  to  the 
bridge  deck.   Resultant  velocities  and  directions  are  derived 
from  these  measurements.   Instruments  are  located  at  several 
stations  along  the  windward  side  of  the  bridge  on  mountings 
which  extend  from  the  edge  of  the  structure  near  roadbed  level. 
Servo  accelerometers  measure  vertical  roadbed  accelerations  with 
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respect  to  time.   The  instruments  are  hard-mounted  on  both  sides 
of  the  bridge  structure  at  roadbed  level.   Acceleration  data  are 
numerically  converted  to  provide  vertical  and  torsional  roadbed 
displacement  data. 

A  block  diagram  of  the  data  recording  system  is  shown  in 
Figure  6.   Analog  signals  from  the  anemometers,  accelerometers 
and  other  sensors  are  continuously  input,  multiplexed  and  con- 
verted to  digital  levels.   These  levels  are  compared  to  preset 
levels  selected  on  the  control  panel,  and  if  greater  for  a 
selected  time  interval,  automatically  activate  data  recording. 
When  the  system  is  active,  continuous  scans  of  all  sensors  are 
input  and  recorded  along  with  coded  time  information  on  digital 
magnetic  tapes.   The  system  has  two  recorders  which  switch  over 
at  the  end  of  a  tape  to  provide  nearly  thirty  hours  of  recording 
capacity.   Tapes  are  standard  nine-track,  800  bpi  NRZI  format 
and  data  are  recorded  in  binary  (uncoded)  form.   Details  of  the 
record  formats  are  included  in  the  data  recording  system  manuals. 

While  overall  system  sampling  occurs  at  500  samples  per 
second,  effective  channel  sample  rates  are  dependant  upon  multi- 
plexer channel  assignments.   At  the  Newport  site,  anemometer 
data  were  sampled  at  14.3  Hz  and  accelerometer  data  were  sampled 
at  7.1  Hz.   At  the  Sitka  site,  both  anemometer  and  accelerometer 
data  were  sampled  at  14.3  Hz.   Eight-bit  analog-to-digital  con- 
version and  signal  conditioning  provided  accelerometer  resolu- 
tion of  .0389  g  at  Newport,  and  .0031  g  at  Sitka.   Anemometer 
resolution  at  each  site  was  .934  mph  (.418  m/s) .   Further  de- 
tails of  the  data  acquisition  system  are  contained  in  other 
publications  (6,7,8,10,12). 

1.2    DATA  REDUCTION 

Data  recorded  on  digital  magnetic  tapes  are  subsequently 
processed  using  FHWA  computer  facilities.   Table  1  lists  many 
of  the  production  oriented  computer  programs  used  to  reduce  the 
data.   Production  software  is  distinguished  from  developmental 
software  in  that  the  former  category  includes  programs  which 
require  little  set-up,  are  seldom  modified  and  which  are  typi- 
cally used  in  processing  project  data.   Developmental  programs 
are  used  for  special  data  reduction  requirements,  and  are  rela- 
tively complex  and  inefficient.   If  a  continuing  need  exists 
for  a  developmental  program,  the  routine  is  revised  to  facili-  , 
tate  production  usage. 

The  data  presented  in  this  report  were  processed  using  the 
spectral  analysis  programs.  These  programs  and  their  output  are 
described  in  detail  in  Sections  2  and  3.   The  output  from  two 
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Table  1.   Significant  project  computer  programs. 


Name 

Purpose 

ACCPLT 

Generates  side-by- s-ide  plots  of  data  from  all 
accelerometers  for  selected  data  records.  . 

ACCSPC 

Calculates  and  plots  power  spectral  density 
for  selected  accelerometer  data. 

CALFIL 

Updates  and  maintains  instrumentation  system 
channel  calibration  file. 

COPY 

Copies  selected  records  from  data  base  to  work 
tape. 

INTERP 

A  subroutine  to  generate  missing  data  points 
through  interpolation. 

MODE  3D 

Plots  three-dimensional  views  of  geometrical 
bridge  configuration  at  sequential  instants  of 
time. 

NEWPRT 

A  subroutine  to  read  Newport  data  tapes,  and 
demultiplex  and  calibrate  sensor  data. 

RBMAX 

Generates  table  of  record-by-record  accelero- 
meter maxima,  and  average  wind  directions  and 
velocities. 

SITKA 

A  subroutine  to  read  Sitka  data  tapes,  and  de- 
multiplex and  calibrate  sensor  data. 

TAP SUM 

Generates  index  and  10-minute  statistical  sum- 
maries for  all  data  recorded  on  field  data  tapes. 

UVCAL 

A  subroutine  to  calculate  anemometer  propeller 
response  correction  factors. 

WINDOW 

2 
A  subroutine  to  apply  10%  cosine   taper  to  data 

samples. 

WNDSPC 

Calculates  and  plots  power  spectral  density  for 
selected  anemometer  data. 

other  programs  provide  unique  insight  to  the  project  from  the 
aspects  of  data  collection  and  data  presentation.  These  pro- 
grams are  briefly  described  in  the  following  paragraphs. 

The  data  acquisition  system  is  activated  automatically 
when  wind  velocities  or  bridge  motions  exceed  certain  thres- 
holds.  Likewise,  the  system  stops  recording  when  the  data  falls 
below  these  levels.   The  obvious  question  then  arises  "What 
data  were  recorded?"   The  tape  summary  programs,  one  for  Newport 
tapes  and  one  for  Sitka  tapes,  were  designed  to  process  entire 
tapes  and  to  summarize  the  events,  conditions  and  data  recorded. 
Figure  7  is  an  example  of  a  statistical  data  summary  output  by 
these  programs.   One  such  summary  is  generated  for  every  ten 
minutes  of  data  recorded  within  a  recording  event.   The  start, 
stop  and  duration  times  for  each  summary  are  indicated  and  the 
minimum,  maximum,  mean  and  standard  deviation  of  all  samples 
for  each  instrument  are  listed.   With  this  output,  it  is  possi- 
ble to  determine  when  the  system  was  recording,  what  activated 
the  system,  the  worst  case  and  general  wind  and  bridge  motion 
activities,  the  similarity  of  activities  between  instruments, 
and  the  integrity  of  the  instruments  and  data  recording  system. 
Foremost,  however,  the  output  permits  the  isolation  of  data  of 
particular  interest,  which  is  a  unique  problem  when  voluminous 
data  is  automatically  recorded  at  remote,  unattended  locations. 

M0DE3D  is  a  computer  program  which  presents  data  of  inter- 
est in  a  form  strictly  for  visual  interpretation.   The  output 
from  this  program  can  be  displayed  on  paper  or  photographic 
film.   Figure  8  shows  several  frames  of  output  which  represents 
in  three-dimensional  perspective  the  geometrical  configuration 
of  the  Sitka  bridge  at  several  instants  of  time.   One  frame  is 
produced  for  each  scan  of  the  accelerometers.   When  viewed  se- 
quentially, these  frames  present  an  animated  reconstruction  of 
the  behavior  of  the  bridge.   Output  to  film  is  effected  using  a 
Stromberg  Carlson  Data  Graphics  SC4020  camera  system  at  Johns 
Hopkins  University  Applied  Physics  Laboratory.   Standard  16mm 
film  is  used  which  can  be  prepared  and  projected  at  speeds  to 
simulate  real-time  or  slow  motion  bridge  behavior.   Viewing 
these  films  provides  unique  insight  into  the  composition  of 
complex  bridge  bending  modes. 

2.     WIND  DATA 

The  wind  data  selected  for  analyses  were  recorded  on  four 
separate  occasions  at  Newport  and  Sitka.   The  most  significant 
data  were  recorded  at  Newport  in  August,  1971  during  the  pass- 
age of  Tropical  Storm  DORIA.   Gusts  up  to  80  mph  accompanied 
the  storm  and  mean  winds  of  nearly  60  mph  were  recorded  across 
the  bridge  deck.   Ten-minute  data  samples  at  various  wind  veloc- 
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STOP  129:?):)9:W 
DURATION     0:  9:56 


SENSOR 

MINIMUM 

MAX  IMUH 

MEAN 

STD.DEV. 

Al 

-0.003 

0.016 

0.0053 

0.0018 

A2 

-0.006 

0.006 

0.0000 

0.0008 

A3 

-0.013 

0.019 

0.0043 

0.3317 

A4 

-0.009 

0.016 

0.0018 

0.0017 

A5 

-0.006 

0.016 

0.0035 

0.0012 

A6 

-0.C09 

0.009 

0.0002 

0.0009 

A7 

-0.013 

0.009 

0.0009 

0.3315 

A8 

-0.009 

0.009 

0.0005 

0.0012 

A9 

-0.009 

0.003 

-0.0037 

0.0014 

A10 

-0.006 

0.006 

0.0001 

0.0008 

All      . 

-0.399 

0.402 

-0.0243 

O.t  33  5 

A12 

-0*031 

0.028 

-0.0004 

0.0052 

A13 

-0.022 

0.028 

-0.0024 

0.0035 

1U 

-10.275 

17.748 

1  .9835 

3.3311 

2U 

-5.605 

10.275 

1.547b 

2.24b0 

3U 

-17. 748 

10.275 

1.6  7b  2 

2.bb51 

4U 

-14.012 

10.275 

0.5507 

3.2004 

5U 

-12.U4 

10.275 

0.8b39 

2.853b 

IV 

0.0 

39.233 

17.9534 

b.78b5 

2V 

2  .802 

35.497 

18.2446 

5.2424 

3V 

3.73b 

39.233 

20.3352 

5.4311 

4V 

1.868 

31.760 

15.7992 

t.mo 

5V 

2.802 

36.4  31 

19.9b70 

5.9127 

1  W 

-7 .473 

14.012 

0.0645 

2.3799 

2W 

-b.539 

11.209 

2.0944 

2.H73 

3W 

-3.736 

9.341 

2.0527 

2.1225 

4W 

-10.275 

13.078 

0.6234 

2.4798 

5W 

-3.73b 

14.012 

2.3521 

2.1214 

SEIS 

0.15b 

0.156 

0.15b3 

0.0  3  00 

SPAR 

-1 .484 

0.0 

-0.5b42 

0.18  37 

SPAR 

-2.109 

-0.625 

-1.1971 

0.1989 

TEMP 

5.800 

b.525 

b.2922 

0.1782 

CL:F 

9.922 

9.922 

9.9219 

0.0000 

CL:Z 

0.0 

0.0 

0.0 

0.0000 

CL:M 

5.000 

5.000 

5.0000 

0.3  3  00 

LV* 

-8.594 

-8.359 

-8.4808 

3.)>iS 

PUM 

0.078 

0.078 

0.0781 

0.0000 

IRG 

0.0 

0.0 

0.0 

0.0000 

Figure  7.   Example  tape  summary  program  statistical  summary, 
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Figure  8.   Example  M0DE3D  program  display. 
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ities  were  processed  to  derive  spectra  which  show  the  general 
amplitude  and  frequency  of  wind  gusts. 

Spectra  for  both  horizontal  and  vertical  winds  were  de- 
rived using  time-to-frequency  domain  transform  techniques.   This 
section  discusses  the  methods  used  to  condition  and  transform 
tne  data,  and  presents  in  separate  subsections,  the  average 
horizontal  and  vertical  wind  spectra  for  each  of  the  four  re- 
cording events.   Additional  detailed  spectra  are  included  in 
Appendices  C-F. 

2.1    PROCESSING  APPROACH 

Alternative  processing  methods  were  tested  and  discussed 
in  a  previous  report  (3).   The  processing  approach  presented  in 
tnis  section  was  adopted  in  consonance  with  the  previous  inves- 
tigation.  Two  additional  considerations  in  conditioning  data 
samples  were  investigated  and  the  results  are  discussed  in 
Sections  2.1.1  and  2.1.2.   Section  2.1.3  describes  the  method 
used  to  correlate  project  data  with  theoretical  models  of  wind 
spectra,  and  lists  the  results  of  the  comparisons. 

Figure  9  shows  the  sequence  of  processing  the  recorded 
wind  data.   Ten-minute  data  samples  of  U-,  V-,  and  W-  component 
winds  for  all  sensor  stations  are  first  isolated  and  corrected 
tor  propeller  response  errors.   Data  sets  are  then  decimated 
to  obtain  a  reasonable  number  of  data  points  and  an  acceptable 
sample  rate.   The  previously  recommended  factor  of  four  is  used 
in  this  study  to  obtain  2048  data  points  at  an  effective  sample 
rate  of  approximately  3.6  Hz.   The  mean  and  linear  ramp  trend, 
if  any,  are  then  removed  and  a  10%  cosine2  taper  is  applied 
prior  to  transforming  the  data.   The  IBM  subroutine  RHARM  is 
used  to  calculate  Fourier  series  coefficients  (5).   These  are 
converted  to  amplitude  and  used  in  the  calculation  of  reduced 
power  density,  RPD,  as  follows: 


.■■■.■■■-■■ 
RPD.  = 


frequency,  x  amplitude. 
i l 

'i    data  sample  variance 


As  a  matter  of  convention,  spectral  frequency  is  converted  to 
wave  number,  WN,  by: 

frequency . 

WN.  =  -= i_ 

i   data  sample  mean  horizontal  velocity 

Nineteen  point  data  averaging  is  used  in  the  final  processing 
step  to  reduce  "spiky"  fluctuations  in  the  spectra  and  to  bet- 
ter display  the  spectral  power  density. 
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TIME  DOMAIN  DATA  SAMPLE 


APPLY  PROPELLER  RESPONSE  CORRECTION 


DECIMATE 


t 

DEMEAN/DETREND 


APPLY  DATA  SAMPLE  TAPER 


TRANSFORM  TO  FREQUENCY  DOMAIN 


COMPUTE  REDUCED  POWER  DENSITY 


CONVERT  FROM  FREQUENCY  TO  WAVE  NUMBER 


SMOOTH  SPECTRAL  DATA 


GENERATE  SPECTRAL  PLOT 


Figure  9.   Functional  processing  sequence  for  wind  data 
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Figure  10.   Example  plot  format  of  wind  spectrum. 


To  simplify  interpretation,  all  wind  spectra  are  displayed 
in  a  common  format.   Figure  10  is  an  example  of  the  plot  format. 
The  abscissa  is  a  three  decade  logarithmic  scale  of  wind  wave 
number  expressed  in  cycles  per  meter. .  The  linear  vertical  scale 
is  reduced  power  density  expressed  in  non-dimensional  units. 
Individual  channel  plots  have  an  inset  listing  the  mean  and  var- 
iance of  the  data  sample  in  miles  per  hour  and  meters  per  second. 
These  plots  also  include  time-traces  of  the  tapered  data  samples 
used  to  derive  the  adjoining  spectra.   The  abscissa  on  these 
traces  is  time  in  4.48  sec  steps,  and  the  ordinate  is  wind  veloc- 
ity in  m/s  from  the  mean. 

2.1.1  Interpolation  of  Discontinuous  Data 

The  data  acquisition  system  installed  at  the  Newport  and 
Sitka  sites  from  October  1970  to  August  197  5  used  incremental 
tape  drives  for  data  recording.   These  instruments  transferred 
measurements  to  magnetic  tape  as  they  were  taken.   When  inter- 
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record  gaps  were  recorded,  however,  measurements  could  not  be 
transferred  and  were  lost.   To  allow  sufficient  time  for  re- 
cording interrecord  gaps,  two  complete  multiplexor  scans  were 
skipped  at  the  end  of  each  record.   This  resulted  in  data  re- 
cording discontinuities  of  .28  sec  after  every  record  (4.48  sec) 
Figure  11  illustrates  the  effect  of  these  discontinuities.   In 
September,  197  5,  the  incremental  tape  drives  were  replaced  with 
fully  buffered,  synchronous  tape  recorders  which  eliminate  data 
discontinuities. 


RECORD  N 


INTER- RECORD 
GAP 


RECORD  N  +  1 


DATA 
POINTS 


SAMPLE  INTERVAL 
.07  sec 


INTERPOLATION  APPLIED 
TO  SHADED  AREA  TO 
SIMULATE  UNRECORDED 
DATA. 


Figure  11.   Recording  system  data  discontinuities. 


To  process  previously  recorded  data,  interpolation  is  used 
to  replace  lost  data  measurements.   For  Newport  recordings,  four 
anemometer  values  and  two  accelerometer  values  are  interpolated. 
For  early  Sitka  recordings,  four  values  are  interpolated  for 
both  anemometers  and  accelerometers.   Aitken's  method  of  inter- 
polation (1)  is  used  to  construct  a  polynomial  of  n-1  degrees 
through  n  given  data  points  which  symetrically  bound  m  missing 
data  points.   The  algorithm  uses  iterative  determinant  solutions 
of  a  polynomial  to  calculate  individual  data  points.   Each  re- 
sulting interpolated  value  has  positional  significance  equal  to 
all  other  points  in  the  data  stream. 

To  test  the  accuracy  of  the  technique,  two  data  sets  were 
generated  to  simulate  real  data.   One  set  approximates  bridge 
motion  data  with  a  .73  Hz  sine  wave.   The  other  data  is  a  set  of 
uniformly  distributed  random  numbers  with  a  mean  of  33.5  and  a 
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standard  deviation  of  5.6,  which  simulates  wind  velocity  data. 
Multiple  data  records  were  created  within  each  set  and  four 
interpolated  data  points  were  calculated  to  correspond  with  four 
data  points  between  adjacent  records.   The  simulated  data  points 
from  the  test  data  and  the  interpolated  points  were  then  separ- 
ately averaged  and  compared.   These  statistics  are  listed  in 
Table  2  and  show  the  relative  accuracy  of  the  interpolation 
technique  as  the  number  of  given  data  points  is  varied. 

Since  interpolated  values  are  determined  by  an  approxima- 
tion of  the  behavior  of  surrounding  points,  accuracy  is  depen- 
dant upon  the  repeatability  and  continuity  of  the  given  data. 
As  expected,  the  errors  associated  with  the  sine  wave  test  data 
are  small  compared  to  the  random  data.   Based  on  these  tests, 
interpolation  using  six  given  data  points  (three  on  each  side  of 
the  data  discontinuity)  provides  maximum  accuracy  for  both  data 
types. 

Table  2.   Expected  errors  for  various  combinations 
of  interpolated  data  points. 


a)   Sine  wave  test  data. 


No. 
Points 

Average 
Test  Value 

Average 
Interp  Value 

Average 
Difference 

Error 

4 

6 

8 

10 

.218140 
.218140 
.218140 
.218140 

.215128 
.217952 
.218133 
.218140 

.003012 

.000188 

.000007 

0.0 

1.3% 
.09% 
.0032% 
0% 

b)   Random 

test  data. 

No. 
Points 

Average 
Test  Value 

Average 
Interp  Value 

Average 
Difference 

Error 

4 

6 

8 

10 

33.786865 
33.786865 
33.786865 
33.786865 

33.559189 
33.639175 
33.156097 
30.647079 

.227676 

.147690 

.630768 

3.139786 

.67% 
.43% 
1.8% 
9.3% 

2.1.2  Correction  for  Anemometer  Response  Characteristics 

An  anemometer  signal  results  from  the  wind  driving  a  poly- 
styrene propeller  which  is  coupled  to  a  miniature  d.c.  voltage 
generator.   The  voltage  is  proportional  to  propeller  speed  and 
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its  polarity  indicates  the  direction  of  propeller  rotation. 
The  propeller  responds  only  to  that  component  of  the  wind  which 
is  parallel  to  its  axis  of  rotation.   Propeller  response  to 
angular  winds  is  a  function  of  the  cosine  of  the  wind  angle. 

Figure  12  compares  the  theoretical  and  the  manufacturer's 
measured  response  data  for  the  four-bladed  propeller  used  in  the 
data  acquisition  system.   By  correcting  for  the  differences  be- 
tween the  calibration  data  and  the  cosine  function,  correction 
factors  can  be  derived  based  on  the  following  relationships: 

9  =  arctan  - 
u 


U  = 


u  cos 


V  =  v  cos  (90-6) 
y90-6 


/U2  +  v2 


where: 


9  is  the  wind  angle  relative  to  the  U-component  axis 
u  and  v  are  recorded  orthogonal  component  velocities 
U,  V,  and  R  are  corrected  component  and  resultant  veloci- 
ties, and 
y  is  the  calculated  response  for  a  given  wind  angle. 

An  eight  degree  polynomial  was  fitted  to  data  points  con- 
sisting of  propeller  responses  and  corresponding  wind  angles 
between  2°  and  88°.   This  limited  data  set  was  used  since  the 
propeller  response  is  symetrical  for  all  positive  direction 
winds  and  since  negative  direction  winds  are  not  considered  in 
project  analyses.   Wind  angles  less  then  2   give  100%  response 
and  wind  angles  greater  than  88   result  in  propeller  stall. 

The  form  of  the  polynomial  and  the  derived  coefficients 
are: 


y  =  a.  +  a,x  +  a_x  +  a^x   ...  a-x 


where : 


y  is  percent  response, 

x  is  wind  angle  in  radians, 
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and  the  following  coefficients  are  given: 


aQ  =  99.70563 
ax  =  5.49133 
a2  =  64.69608 
a.-  =  69.5257  4 


a.  =  85.19067 
4 

ac  =  25.50650 

D 

ac   =  56.69714 
6 

a?  =  17.33202 


Table  3  lists  the  results  of  propeller  response  correction  for 
various  wind  angles.   All  figures  are  percent  of  response  to 
resultant  wind  velocity. 

Table  3.   Corrected  anemometer  responses  and 
remaining  errors  at  selected  wind  angles. 


Wind 

Ideal 

Actual 

Corrected 

Remaining 

Angle 

Response 

Response  . 

Response 

Error 

(degrees) 

(%) 

(%)  (a) 

(%)  (b) 

(%)  (a-b) 

0 

100.0 

100.0 

100.0 

0.0 

10 

98.5 

98.0 

98.1 

0.4 

20 

94.0 

92.0 

93.9 

0.1 

30 

86.6 

82.0 

87.4 

-0.8 

40 

76.6 

68.0 

77.0 

-0.4 

50 

64.3 

51.0 

61.9 

2.4 

60 

50.0 

40.0 

52.1 

-2.0 

70 

34.2 

24.0 

34.0 

-0.2 

80 

17.4 

10.0 

16.8 

-0.6 

90 

0.0 

0.0 

0.0 

0.0 

As  an  example,  given  a  50  mph  wind  at  an  angle  of  30   from  the 
U-axis,  the  ideal  U-component  velocity  should  be  86.6%  x  50  = 
43.3  mph.   However,  the  measured  velocity  after  correction  will 
be  87.4%  x  50  =  43.7  mph,  resulting  in  an  error  of  0.8%  x  50  = 
.4  mph.   The  mean  error  over  all  angles  is  0.05%  with  a  standard 
deviation  of  1.05%. 

2.1.3  Correlation  of  Data  with  Theoretical  Models 

The  plots  in  Sections  2 . 2  and  2 . 5  show  the  combined  spec- 
tral results  for  all  transducers  for  each  particular  data  sample. 
Average  horizontal  and  vertical  spectra  are  included  for  each 
selected  nominal  wind  velocity.   An  approach  to  superimpose  the 
curves  of  theoretical  models  on  these  spectra  was  described  in 
a  paper  coauthor ed  by  the  FHWA  contract  manager,  Mr.  R.H.  Gade 
(2) .   An  extract  from  this  paper  is  included  below: 
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For  horizontal  spectra  a  comparative  plot  is  based  on  the 
following  Harris  -  Von  Karman  formula:   (4) 


nSu(n) 
kU2  (10) 


4N 


(2+Nz) 


2^5/6 


where  k  is  the  surface  drag  coefficient,  and  U(10)  is  the  mean 
velocity,  both  referred  to  10-meter  height,  and 


N  = 


nL 


U(10) 


,  L  being  a  peak  wavelength  "constant". 


This  formula  is  converted  for  comparative  purposes  to  the  de- 
sired form  as  follows: 


nSyCn) 


U 


4K^ir  (z) 


2  InM  m  62.7 


U 


Zo 


z0 


Hz) 


In 


62.7 

Zo 


5/6 


2  + 


Mz) 


In 


10 

Zo 


where  K,  Karman 's  constant,  is  taken  as  0.4,  U(z)  is  averaged 
mean  wind  velocity  at  bridge  height  and  1/A(z)  =  n/U(z);  z0, 
the  surface  roughness  length,  is  related  to  k  by  the  formula: 


K  = 


K 


(*  ^ 


The  peaks  of  the  sample  and  reference  spectra  were  arbi- 
trarily equated  to  obtain  the  values  of  k  and  L.   The  resulting 
relations  are: 


1.812  kU  (10)  =  [n  Su(n)l 

/I  m(10/z°) 


max 


L  = 


A(z) 


In 


62.7 

zo 
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For  vertical  spectra  an  empirical  spectrum  is  derived 
from  the  following  Panofsky  formula:   (9) 


nSW(n)        6£ 


kU2 (10)   (1  +  4f) 


where  f  =  nz/U,  z  being  the  height  of  the  transducer,  and  U 
being  the  mean  horizontal  velocity  at  that  height. 

This  formula  is  converted  for  comparative  purposes  to 
the  desired  form  as  follows: 


nSW(n)         K2  U2(z) (l/A(z))  6z 


W 


0w(^^)2(1  +  4z/X<Z) 


The  surface  drag  coefficient  is  similarly  obtained  by 

0.375  kU2  (10)  =  [nST7(n)l 

L  W   J  max 

and  zo  as  above. 

Peak  values  from  sample  spectra  were  identified  manually 
and  a  computer  program  was  used  to  generate  smooth  curves  and 
calculate  implied  model  parameters  based  upon  these  values. 
The  peak  values,  model  parameters  and  sample  identification  in- 
formation are  listed  in  Table  4. 
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2.2    NEWPORT  -  AUGUST  1971 

This  section  presents  data  recorded  at  the  Newport  bridge 
site  on  August  28,  1971.   The  specific  start  and  stop  times  asso- 
ciated with  the  data  are  unknown  due  to  recorded  time  inaccura- 
cies.  The  data  are  from  Tropical  Storm  DORIA,  which  passed  near 
the  bridge,  with  mean  wind  velocities  to  60  mph  and  gusts  to  80 
mph.   Data  from  five  wind  sensors  were  processed  at  mean  wind 
velocities  of  20,  35,  50,  and  60  mph. 


Figure  13  compares  the  average  vertical  wind  angles  of 
attack  at  the  mean  U-component  wind  velocities  experienced  at 
each  sensor  station.   Figures  14  through  17  show  five-sensor, 
ensemble  average  wind  spectra  at  the  four  selected  wind  veloci- 
ties.  Each  figure  includes  separate  plots  of  the  U-component 
spectrum,  the  resultant  horizontal  spectrum,  the  vertical  spec- 
trum adjusted  to  the  U-component  velocity,  and  the  vertical 
spectrum  adjusted  to  the  resultant  horizontal  velocity.   Indivi- 
dual channel  spectra  for  this  recording  event  are  included  in 
Appendix  C. 
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Figure  13.   Angle  of  attack  vs.  wind  velocity  at  individual 
sensor  locations  (Newport,  August  1971) . 
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k   i  i  mag    *   i  i  romp — i   i  i  romp  T~n  »n» 

Wave  No . 
a)    U-component  -    5U,    6U,    7U,    8U,    9U. 


i     4  4  4  4444lf» — 4     4  4  4  4UW — 4     4   4  lUW 

Wave  No . 


4  4  4  4444 


b)  Resultant  horizontal  -  5R,  6R,  7R,  8R,  9R. 

Figure  14.   Five-channel  average  horizontal  and  vertical 
wind  spectra  for  20  mph  wind  (Newport,.  August  1971) . 
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Wave  No. 
c)  Vertical  -  5W/U,  6W/U,  7W/U,  8W/U,  9W/U. 


d)  Vertical  -  5W/R,  6W/R,  7W/R,  8W/R,  9W/R. 
Figure  14.   Five-channel  average  horizontal  and  vertical  wind 
spectra  for  20  mph  wind  (Newport,  August  19-71)  (continued)  . 
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b)  Resultant  horizontal  -  5R,  6R,  7R,  8R,  9R. 
Figure  15.   Five-channel  average  horizontal  and  vertical 
wind  spectra  for  35  mph  wind  (Newport, .  August  1971). 
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Wave  No . 
c)  Vertical  -  5W/U,  6W/U,  7W/U,  8W/U,  9W/U. 


v    *  rrromy    i  i  rrnng    i  rTrettg — i  i  rrnw 

Wave  No. 

d)  Vertical  -  5W/R,  6W/R,  7W/R,  8W/R,  9W/R. 
Figure  15.   Five-channel  average  horizontal  and  vertical  wind 
spectra  for  35  mph  wind  (Newport,  August  19.71)  (continued)  . 
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Wave  No. 
a)     U-component    -    5U,    6.U,    7U,    8.U,    9U. 


Wave  No . 

b)  Resultant  horizontal  -  5R,  6R,  7R,  8R,  9R. 
Figure  16.   Five-channel  average  horizontal  and  vertical 
wind  spectra  for  50  mph  wind  (Newport,.  August  1971) . 
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Wave  No. 
c)  Vertical  -  5W/U,  6W/U,  7W/U,  8W/U,  9W/U. 


mur 


Wave  No. 


rmw 


d)  Vertical  -  5W/R,  6W/R,  7W/R,  8W/R,  9W/R. 
Figure  16.   Five-channel  average  horizontal  and  vertical  wind 
spectra  for  50  mph  wind  (Newport,  August  1971) (continued) . 
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*l? — I     *  4  WWmSEF     I     *  4  imUgi — 1     4  4  U\UW — I     4  4  14*41 


Wave   No. 
a)    U-component   -    5U,    6U,    7U,    8U,    9U. 


Wave  No . 

b)  Resultant  horizontal  -  5R,  6R,  7R,  8R,  9R. 
Figure  17.   Five-channel  average  horizontal  and  vertical 
wind  spectra  for  60  mph  wind  (Newport,  August  1971). 
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Wave  No. 
c)  Vertical  -  5W/U,  6W/U,  7W/U,  8W/U,  9W/U. 


»b*    i   rrrnag    *  i  hmmw    i   r-rnrmg — i   i  hum 

Wave  No. 

d)  Vertical  -  5W/R,  6W/R,  7W/R,  8W/R,  9W/R. 
Figure  17.   Five-channel  average  horizontal  and  vertical  wind 
spectra  for  60  mph  wind  (Newport,  August  19.71)  (continued)  . 
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2.3    NEWPORT  -  JANUARY  1972 

This  section  presents  data  recorded  at  the  Newport  bridge 
site  in  January  1972.   The  specific  date,  start  and  stop  times 
associated  with  the  data  are  unknown  due  to  recorded  .time  inac- 
curacies.  The  data  are  from  a  storm  with  mean  wind  velocities 
to  50  mph  and  gusts  to  70  mph.   Data  from  five  wind  sensors  were 
processed  at  mean  wind  velocities  of  40  and  50  mph. 

Figure  18  compares  the  average  vertical  wind  angles  of 
attack  at  the  mean  U-component  wind  velocities  experienced  at 
each  sensor  station.   Figures  19  and  2  0  show  five-sensor,  en- 
semble average  wind  spectra  at  the  two  selected  wind  velocities. 
Each  figure  includes  separate  plots  of  the  U-component  spectrum, 
the  resultant  horizontal  spectrum,  the  vertical  spectrum  adjust- 
ed to  the  U-component  velocity,  and  the  vertical  spectrum  ad- 
justed to  the  resultant  horizontal  velocity.   Individual  channel 
spectra  for  this  recording  event  are  included  in  Appendix  D. 
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Figure  18.   Angle  of  attack  vs.  wind  velocity  at  individual 
sensor  locations  (Newport,  January  1972) . 
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b)  Resultant  horizontal  -  2R,  3R,  5R,  6R,  7R. 
Figure  19.   Five-channel  average  horizontal  and  vertical 
wind  spectra  for  40  mph  wind  (Newport,  January  1972) . 
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Wave  No . 
c)  Vertical  -  2W/U,  3W/U,  5W/U,  6W/U,  7W/U. 


Wave  No . 


ITTT 


d)  Vertical  -  2W/R,  3W/R,  5W/R,  6W/R,  7W/R. 
Figure  19.   Five-channel  average  horizontal  and  vertical  wind 
spectra  for  4  0  mph  wind  (Newport,  January  197  2) (continued) . 
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Wave  No. 
a)    U-component  -   2U,    3U,    5U,    6U,    7U. 


»U»   i  rrrniny   i  rrrnn&   i  rrmm? — i  i  rmw 

Wave  No. 

b)  Resultant  horizontal  -  2R,  3R,  5R,  6R,  7R. 
Figure  20.   Five-channel  average  horizontal  and  vertical 
wind  spectra  for  50  mph  wind  (Newport,  January  1972) . 
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Wave  No . 
c)  Vertical  -  2W/U,  3W/U,  5W/U,  6W/U,  7W/U. 
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I  iuui 


t  4  nmtf 

Wave  No. 
d)  Vertical  -  2W/R,  3W/R,  5W/R,  6W/R,  7W/R. 

Figure  20.   Five-channel  average  horizontal  and  vertical  wind 
spectra  for  5  0  mph  wind  (Newport,  January  1972) (continued) . 
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2.4    SITKA  -  NOVEMBER  1975 

This  section  presents  data  recorded  at  the  Sitka  bridge 
site  at  2218  on  November  10,  1975.   The  processed  data  have  a 
mean  wind  velocity  of  30  mph  and  gusts  to  50  mph.   Data  from 
four  wind  sensors  were  processed,   Since  the  wind  direction  was 
nearly  parallel  to  the  bridge,  only  the  resultant  horizontal 
wind  was  used.   U-component  winds  were  not  separately  displayed, 

Figure  21  compares  the  average  vertical  wind  angles  of 
attack  at  the  mean  resultant  horizontal  wind  velocities  experi- 
enced at  each  sensor  station.   Figure  22  shows  four-sensor,  en- 
semble average  wind  spectra  at  the  selected  wind  velocity.   The 
plot  at  the  top  of  the  figure  is  the  resultant  horizontal  wind 
spectrum  and  the  plot  at  the  bottom  is  the  vertical  wind  spec- 
trum.  Individual  channel  spectra  for  this  recording  event  are 
included  in  Appendix  E. 
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Figure  21.   Angle  of  attack  vs.  wind  velocity  at  individual 
sensor  locations  (Sitka,  November  1975) . 
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Wave  No . 
a)  Resultant  horizontal  -  2R,  3R,  4R,  5R. 


Wave  No . 
b)  Vertical  -  2W/R,  3W/R,  4W/R,  5W/R. 

Figure  22.  .  Four-channel  average  horizontal  and  vertical 
wind  spectra  for  3  0  mph  wind  (Sitka,  November  197  5) . 
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2.5    SITKA  -  MAY  1976 

This  section  presents  data  recorded  at  the  Sitka  bridge 
site  at  0536  on  May  9,  1976.   The  processed  data  have  a  mean 
wind  velocity  of  25  mph  and  gusts  to  35  mph.   Data  from  four 
wind  sensors  were  processed.   Since  the  wind  direction  was  near- 
ly parallel  to  the  bridge,  only  the  resultant  horizontal  wind 
was  used.   U-component  winds  were  not  separately  displayed. 

Figure  23  compares  the  average  vertical  wind  angles  of 
attack  at  the  mean  resultant  horizontal  wind  velocities  experi- 
enced at  each  sensor  station.   Figure  24  shows  four-sensor  en- 
semble average  wind  spectra  at  the  selected  wind  velocity.   The 
plot  at  the  top  of  the  figure  is  the  resultant  horizontal  wind 
spectrum,  and  the  plot  at  the  bottom  is  the  vertical  wind  spec- 
trum.  Individual  channel  spectra  for  this  recording  event  are 
included  in  Appendix  F. 
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Figure  23.   Angle  of  attack  vs.  wind  velocity  at  individual 
sensor  locations  (Sitka,  May  1976) . 
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Wave  No . 
a)    Horizontal   -    1R,    2R,    3R,    5R. 


nir 


Wave  No. 

b)  Vertical  -  1W/R,  2W/R,  3W/R,  5W/R. 

Figure  24.   Four-channel  average  horizontal  and  vertical 
wind  spectra  for  25  mph  wind  (Sitka,  May  1976) . 
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3.     BRIDGE  MOTION  DATA 

The  bridge  motion  data  selected  for  analyses  were  recorded 
at  Sitka  on  a  single  occasion.   No  other  bridge  motion  data  are 
available  due  to  either  instrumentation  problems  or  to  lack  of 
significant  bridge  motion  activity. 

In  this  section,  acceleration  data  are  presented  in  the 
form  of  time  history  traces  and  spectral  plots.   The  methods  used 
to  process  the  data  and  the  results  of  processing  are  discussed. 
Detailed  spectral  plots  are  included  in  Appendix  G. 

3.1    PROCESSING  APPROACH 

Figure  25  indicates  the  sequence  of  steps  used  to  process 
accelerometer  data.   This  sequence  is  the  same  as  that  used  to 
process  anemometer  data  except  for  sample  preparation  and  the 
conversion  of  terms  for  data  presentation.   Continuous  512-point 
data  samples  are  used  to  obtain  maximum  discrete  frequency  in- 
formation in  distinct  bridge  motion  events.   Samples  are  demean- 
ed and  detrended  to  eliminate  bias  and  linear  drift  errors,  and 
are  then  tapered  and  transformed  into  the  frequency  domain. 
Conventional  power  density  values  are  derived  from  Fourier  ser- 
ies coefficients  and  the  data  are  plotted  with  respect  to  fre- 
quency. 

An  example  plot  is  shown  in  Figure  26.   The  vertical  axis 
is  power  density  in  units  of  g^/sec.   The  range  of  this  axis 
varies  from  plot  to  plot  to  show  the  relative  amplitude  of  fre- 
quency components  within  the  data  sample.   The  horizontal  axis 
is  expressed  in  Hertz.   A  time  history  of  the  tapered  data  sam- 
ple is  shown  in  the  inset  plot  accompanying  individual  channel 
spectra.   The  ordinate  is  expressed  in  g's  and  the  abscissa  shows 
relative  time  where  each  mark  represents  4.48  seconds  (one  data 
record) . 

To  distinguish  between  frequency  components  representing 
vertical  and  torsional  bending  modes,  acceleration  data  from 
instruments  on  opposite  sides  of  the  bridge  are  arithmetically 
summed  and  differenced,  and  spectra  are  produced  to  indicate 
whether  the  outputs  of  the  accelerometer  pairs  are  in-phase.   The 
hypothesis  is  that  in-phase  data  implies  vertical  bridge  motion 
and  out-of -phase  data  implies  torsional  bridge  motion.   Since 
frequency  components  in  spectra  of  in-phase  paired  accelerometer 
data  will  be  of  greater  amplitude  than  frequency  components  in 
spectra  of  individual  accelerometer  data,  those  frequencies  with 
greater  amplitudes  must  represent  vertical  bridge  motion.   Simi- 
larly, frequency  components  related  to  torsional  bridge  motion 
will  be  of  lower  amplitudes  in  spectra  of  paired  accelerometer 
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TIME  DOMAIN  DATA  SAMPLE 


DEMEAN/DETREND 


APPLY  DATA  SAMPLE  TAPER 


1 ' 
TRANSFORM  TO  FREQUENCY  DOMAIN 


COMPUTE  POWER  DENSITY 


GENERATE  SPECTRAL  PLOT 


Figure  25.   Functional  processing  sequence 
for  acceleration  data. 
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data,  compared  to  spectra  of  individual  accelerometer  data. 
The  latter  components  are  shown  by  producing  spectra  of  the 
arithmetic  difference  data  between  accelerometer  pairs.   This 
causes  torsional  frequency  components  to  increase  rather  than 
decrease  in  amplitude  and  simplifies  comparisons  of  spectra. 
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Figure  26.   Example  plot  format  of  acceleration  spectrum. 


3.1.1  Calculation  of  Damping  Factors 

Attempts  to  calculate  bridge  structure  damping  factors 
using  direct  measurements  of  the  decaying  time  history  ampli- 
tudes proved  inadequate.   The  calculated  values  were  simply  too 
low  to  be  believable.   An  alternate  method,  however,  provides 
factors  which  are  both  believable  and  consistent.  This  approach 
involves  the  measurement  of  the  bandwidth,  Aoj,  at  the  half -power 
points  of  a  discrete  frequency  component,  w  .   These  measure- 
ments are  used  to  calculate  the  system  damping  factor,  £ ,  by: 
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Acceleration  spectral  data  values  are  plotted  and  measure- 
ments are  taken  manually.   Figure  27  is  a  typical  plot  of  the  • 
discrete  data  points  used  to  define  a  frequency  component  at 
approximately  2.65  Hz.   The  bandwidth  of  this  component,  .06  Hz, 
is  measured  at  the  indicated  half-amplitude  point  and  results  in 
a  damping  factor  of  .011. 
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Figure  27.   Damping  factor  calculations  for  a  typical 
acceleration  spectrum  frequency  component. 


Application  of  this  method  results  in  calculated  damping 
factors  between  .006  and  .011  for  torsional  vibrations  and  .009 
to  .013  for  vertical  vibrations.   These  preliminary  results  are 
based  on  limited  acceleration  data;  further  investigations  will 
be  based  upon  derived  displacement  data. 


3.2 


SITKA  -  DECEMBER  197  5 


Significant  bridge  motion  activity  occurred  during  a  re- 
cording event  on  December  6/  1975.   The  recording  system  was 
activated  by  winds  which  reached  a  peak  velocity  of  38  mph.   The 
system  was  active  from  1426  to  1432.   Several  other  recording 
events  had  occurred  during  the  day  with  similar  winds  but  lower 
amplitude  bridge  vibrations.   The  cause  of  the  bridge  vibration's 
has  not  been  determined.   Figure  28  shows  a  time  history  of  data 
from  all  accelerometers  during  the  recording  event.   The  ordi- 
nate for  the  thirteen  accelerometer  traces  is  amplitude  in  g's 
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a)   First  bridge  motion  event. 

Figure  28.   Time  history  of  acceleration  data 
(Sitka,  December  1975) . 
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b)   Second  bridge  motion  event. 

Figure  28.   Time  history  of  acceleration  data 
(Sitka,  December  1975) (continued) . 
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c)   Third  bridge  motion  event. 

Figure  28.   Time  history  of  acceleration  data 
(Sitka,  December  1975) (continued) . 
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and  the  abscissa  is  time.   The  marks  at  the  top  of  the  plots  in- 
dicate records  and  the  marks  at  the  bottom  of  the  plots  indicate 
seconds.   Three  separate  vibration  events  were  defined  within  ' 
distinct  eight  record  (36  seconds)  intervals.   Data  from  each  of 
these  events  were  processed  separately  because  each  shows  some- 
what different  vibrational  activity. 

Frequency  data  for  the  events  are  illustrated  in  Figures 
29  through  31.   Each  row  of  four  plots  shows  simplified  spectra 
for  the  two  accelerometers  on  opposite  sides  of  the  bridge  at  a 
particular  station;  the  difference  or  out-of -phase  spectrum 
which  o-orresponds  to  torsional  bending  modes;  and  the  sum  or 
in-phase  spectrum  which  corresponds  to  vertical  bending  modes. 
The  plots  show  the  predominant  frequency  components  for  each 
spectrum.   The  detailed  spectra  for  the  accelerometers  are  in- 
cluded in  Appendix  G.   Note  that  the  ordinate  range  which  is 
expressed  in  g^/sec  x  10~6  is  constant  for  each  row  but  varies 
from  row  to  row  and  from  event  to  event.   Variation  in  the 
scale  permits  display  of  predominant  frequency  components  while 
retaining  a  measure  of  absolute' amplitudes. 

In  spectra  for  the  three  events,  some  frequency  components 
increase  and  some  decrease  in  both  relative  and  absolute  ampli- 
tude.  Also,  in  comparing  events,  there  is  some  variability  in 
spectral  frequency  content  and  implied  bending  modes.   However, 
based  upon  the  hypothesis  of  correlation  between  data  from  the 
two  sides  of  the  bridge,  the  predominant  roadbed  torsional  mode 
frequencies  are  2.65  and  5.7  Hz.   The  predominant  roadbed  ver- 
tical mode  frequency  is  2.65  Hz.   Several  strong  in-phase  tower 
acceleration  frequencies  were  present  near  3.5  and  4.2  Hz  and  a 
strong  out-of-phase  frequency  occurred  at  5.7  Hz.   Table  5  pre- 
sents a  more  complete  list  of  predominant  frequencies.   Note 
that  all  frequency  information  is  based  on  acceleration  data. 
Utilization  of  displacement  data  may  yield  slightly  different 
results. 

4.  .    SUMMARY  AND  PLANS 

During  this  reporting  period,  three  trips  were  made  to 
Sitka,  Alaska,  to  calibrate  the  instrumentation.   During  one 
trip,  the  incremental  drive  tape  recorders  were  replaced  with 
buffered,  synchronous  drive  recorders.   The  performance  of  the 
instrumentation  of  the  new  recorders  has  been  excellent  and 
subsequent  data  recordings  have  been  virtually  error  free. 
While  substantial  data  have  been  recorded,  most  winds  have  been 
from  the  southeast  at  moderate  speeds.   Since  the  bridge  has  an 
east-west  alignment,  large  wind  angles  have  been  experienced 
rather  than  the  anticipated  zero-angle,  across-deck  winds. 
Only  one  bridge  motion  event  was  recorded  and  it  may  have  been 
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Figure  29.  Simplified  individual,  sum  and  difference  spectra 
for  first  bridge  motion  event. 
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Figure  31.   Simplified  individual,  sum  and  difference  spectra  for 
third  bridge  motion  event. 
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Table  5.   Summary  of  predominant  spectral 
components  during  bridge  motion  events. 


Event 

Roadbed 

Towers 

Vertical 

Span 

Torsion 

Span 

In-Phase 

Out-Phase 

1 

5.3 

Mid 

5.7 
6.9 

Qtr 
Mid 

3.9 
5.7 

5.4 
5.7 
6.9 

2 

2.5 

End 

2.6 

Qtr 

3.4 
3.6 
3.8 

2.6 

3 

2.6 

End 

2.6 
3.9 
3.9 

Qtr 
Mid 
End 

3.6 
3.9 

All  frequency  values  expressed  in  Hertz. 


caused  by  vehicular  traffic.   Current  plans  are  to  remove  the 
system  from  the  Sitka  Harbor  bridge  in  the  Spring,  1977,  and  re- 
turn it  for  overhaul  at  the  Federal  Highway  Research  Station  in 
McEean,  Virginia.   Other  bridge  sites  are  under  evaluation  as 
candidates  for  future  instrumentation  installations. 

All  data  recorded  at  the  Newport  site,  and  data  recorded 
at  the  Sitka  site  since  the  installation  of  the  new  recorders, 
have  been  cataloged.   A  summary  of  the  data  base  is  included  in 
Appendix  B.   This  report  has  presented  selected  data  from  the 
data  base,  indicative  of  wind  and  bridge  motion  activity  at  the 
two  instrumented,  long-span  bridge  sites.   Frequency  spectra  of 
wind  data  have  been  generated  to  show  the  statistical  charac- 
teristics of  winds  at  various  velocities  at  the  two  sites. 
Frequency  spectra  of  bridge  motion  data  show  the  predominant 
vibratory  frequencies  and  bending  modes  of  the  Sitka  bridge. 
Future  efforts  will  involve  the  calculation  of  bridge  displace- 
ment data  from  acceleration  data.   Corresponding  frequency 
spectra  will  be  compared.   Future  motion  pictures  of  bridge 
motion  will  also  be  based  on  displacement  data.   Should  a  seis-' 
mic  event  or  significant  winds  initiate  bridge  activity  at 
Sitka,  all  efforts  will  be  directed  towards  the  detailed  anal- 
ysis and  correlation  of  bridge  motions  and  the  initiating  force. 
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APPENDIX  A 
CHRONOLOGY  OF  SIGNIFICANT  PROJECT  EVENTS 

06/30/69   Contract  awarded  to  Melpar  (American  Standard, 
Inc.)  to  design, , fabricate,  and  install  the 
instrumentation  system. 

09/20/69   Survey  of  Newport,  Rhode  Island,  bridge  site. 

10/26/70   System  installed  on  Newport  Bridge. 

12/28/70   Calibration  and  maintenance  trip. 

01/28/71  Accelerometers  removed  and  returned  to 
manufacturer  for  repair. 

05/19/71   Sensor  cables  removed  from  bridge  due  to 
unscheduled  painting  of  bridge. 

07/01/71   Contract  awarded  to  Systems  Technology 

Associates,  Inc.  to  process  data  and  ser- 
vice instrumentation  system. 

08/16/71   Reinstall  instrumentation  on  west  half  of 
bridge. 

08/28/71  Tropical  Storm  DORIA. 

08/31/71  Trip  to  Newport  bridge  to  retrieve  data  tapes. 

01/16/72  Reinstall  instrumentation  on  east  half  of  bridge 

02/09/72  Trip  to  bridge  to  retrieve  data  tapes. 

02/24/72   Calibration  and  maintenance  trip.   Shelter 
constructed  for  instrumentation. 

03/20/72   Calibration  and  maintenance  trip. 
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04/24/72   Sensors  disconnected  and  system  shut  down 
for  bridge  repainting. 

06/04/72   System  removed  from  Newport  bridge  site,  and 
returned  to  FHRS  for  refurbishment.  . 

09/30/72   System  installed  on  Sitka,  Alaska  bridge. 

11/09/72  Alaska  Department  of  Highways  replaced    -%• 
accelerometer  A10. 

02/12/73  Calibration  and  maintenance  trip.   Recorder 
#1  removed  and  returned  to  manufacturer 
for  repair. 

07/09/73  Calibration  and  maintenance  trip.   Recorder 
#1  reinstalled. 

09/18/73  Presentation  on  .project  given  at  annual 

conference  on  Federally  Coordinated  Program 
of  R&D  in  Highway  Transportation,  San 
Francisco. 

10/09/73  Trip  to  bridge  site.   Tower  motion  study. 

04/15/74  Calibration  and  maintenance  trip.   Installed 
line- voltage  and  power-up  monitors. 

08/15/74  Calibration  and  maintenance  trip.   Replaced 
accelerometers  A3,  A4. 

11/ — /74  Automobile  accident  on  bridge  damaged  sensor 
cables. 

12/16/74   Calibration  and  maintenance  trip.   Repaired 
cables  from  automobile  accident.   Replaced 
accelerometers  Al,  A2. 

05/11/75   Calibration  and  maintenance  trip.   Refur- 
bished anemometers,  relocated  accelero- 
graph,  replaced  accelerometer  A9. 
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08/24/75  Calibration  and  maintenance  trip.  Instal- 
led new  tape  recorders.  Replaced  acceler- 
ometers  A2,  A4,  A8,  All.  Replaced  anemom- 
eters at  W3,  W4,  W5. 

09/19/75   Presentation  on  project  given  at  annual  con- 
ference on  Federally  Coordinated  Program  of 
R&D  in  Highway  Transportation,  Minneapolis.. 

10/29/75   Calibration  and  maintenance  trip.   Repaired 
tape  recorder  #1  and  faulty  sensor  cables. 
Replaced  anemometers  at  W2 . 

03/19/76   Calibration  and  maintenance  trip. 
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APPENDIX  C 
INDIVIDUAL  CHANNEL  WIND  SPECTRA.:   NEWPORT,  AUGUST  1971 
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NAVE  NO 

b)   Resultant  horizontal  -  5R. 

Figure  32.   Individual  channel  horizontal  and  vertical 
wind  spectra  for  20  mph  wind  (Newport,  August  1971)  . 

67 


MPfl  M/5 

MEflN=U.72     2.11 
VRRI*«1.32     0.86 


MRVE  M 


c)    Vertical   -    5W/U. 
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MEAN- 22. 62  10.11 
VRRI=i6.96  S.39 


WWE  Nfl 


d)  U~component  -  6U. 
Figure  32.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  20  mph  wind  (Newport,  August  1971) (continued) . 
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e)    Resultant   horizontal   -    6R. 
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MERN=5.21     2.33 
VfiRI  =  (4.20     0.8<1 


HAVE  Nil 

f)  Vertical  -  6W/U. 
Figure  32.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  20  mph  wind  (Newport,  August  1971) (continued) . 
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g)  U-component  -  7U. 
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mm  m 

b)  Resultant  horizontal  -  7R. 

Figure  32.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  20  mph  wind  (Newport,  August  1971) (continued) . 
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i)  Vertical  -  7W/U. 
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j)  U-component  -  8U. 

Figure  32.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  20  mph  wind  (Newport,  August  1971) (continued) . 
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k)  Resultant  horizontal  -  8R. 
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MERN-U.89  2.18 
VRRI-U.33  0.87 


iwve  NO 
1)  Vertical  -  8W/U. 

Figure  32.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  20  mph  wind  (Newport,  August  1971) (continued) . 
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m)  U-component  -  9U. 
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NAVI  Nf 

n)  Resultant  horizontal  -  9R. 

Figure  32.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  20  mph  wind  (Newport,  August  1971) (continued) . 
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o)  Vertical  -  9W/U. 

Figure  32.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  20  mph  wind  (Newport,  August  1971) (continued) . 
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MEAN- HO. 96  18.31 
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a)  U-component  -  5U. 
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VRRI-33.98  6.79 
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b)  Resultant  horizontal  -  5R. 
Figure  33.   Individual  channel  horizontal  and  vertical 
wind  spectra  for  35  mph  wind  (Newport,  August  1971) . 
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mm 

MEAN-  U.  60  2.06 
VRRI-6.29  1.26 


c)  Vertical  -  5W/U. 
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»to» — 4    4  4  UUHg* — i    4  4  iUtHg* — i    4  4  4  4444W* — 4    4  4  HU1 

d)  U-component  -  6U. 

Figure  33.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  35  mph  wind  (Newport,  August  1971) (continued) . 
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N/1 
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VRRI-22.V7  V.V9 


e)    Resultant   horizontal    -    6R 


MPH  N/S 

MERN-5.VV    2.V3 
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To" — *   i  i  rnmy    *   i  i  nu^ — *   i  unu^ — i   i  rrnn 

MOVE  NO 

f)  Vertical  -  6W/U. 
Figure  33.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  35  mph  wind  (Newport,  August  1971) (continued) . 
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N/t 

MERM-38.16  17.00 
VflRI-23.6«  «.72 


h)  Resultant  horizontal  -  7R. 
Figure  33.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  35  mph  wind  (Newport,  August  1971) (continued) . 
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N€ftN-5.22    2.33 
VRRI-7.06    1.U1 
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i)  Vertical  -  7W/U. 
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j)  U-component  -  8U. 

Figure  33.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  35  mph  wind  (Newport,  August  1971) (continued) . 
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k)    Resultant   horizontal    -    8R. 
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VRRI-7.58    1.51 


w — i   i  tuuiu*    i   i  i  tmiu" — i   i  iimw — i   t  i  um 

1)  Vertical  -  8W/U. 
Figure  33.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  35  mph  wind  (Newport,  August  1971) (continued) . 
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tut 
HERN- 37. 56  16.79 
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m)  U-component  -  9U. 


NPH  H/S 

HERN* 46. 16  20.64 
VRRI-24.37  4.87 


JU*    *   iiuuiu*    i   itumu*    *   mmw — *   i  iwm 


n)  Resultant  horizontal  -  9R. 

Figure  33.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  35  mph  wind  (Newport,  August  1971) (continued) . 
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2.«« 
0.97 


V   *  rmug   i  rnrrgg    i  tm 


<   *  Uiiir 


o)  Vertical  -  9W/U. 
Figure  33.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  35  mph  wind  (Newport,  August  1971)  (continued) . 
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HAVE  N9 

a)  U-component  -  5U. 


N/S 

HEflN>5U.OO  2U.1U 
VfiRI-34.70  6.93 


b)  Resultant  horizontal  -  5R. 

Figure  34.   Individual  channel  horizontal,  and  vertical 
wind  spectra  for  50  mph  wind  (Newport,  August  1971) . 
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c)  Vertical  -  5W/U. 
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VflRI-35.79  7.15 


»n* — i   h  i  nuw — i   i  i  uuhj* — *   t  i  iitnk< — i   i  *  row 

MfWC  M 

d)  U-component  -  6U. 
Figure  34.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  50  mph  wind  (Newport,  August  1971) (continued) . 
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N/S 

MEAN- 49. 96  22.31 
VRRI-32.37  6.147 


e)    Resultant   horizontal    -    6R. 


NPN  M/S 

HERN- 5. 59     2.50 
VRRI-7.2U     1.U5 


f)  Vertical  -  6W/U. 

Figure  34.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  50  mph  wind  (Newport,  August  1971) (continued) . 
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g)    U-component   -    7U. 
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HkiHikA.^iiiiaj  nr4 


MEAN- 50. 93  22.76 
VRRI-29.91  5.96 
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5*     *     4  4  tIMJk* 


h)  Resultant  horizontal  -  7R. 

Figure  34.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  50  mph  v/ind  (Newport,  August  1971) (continued) . 
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i)  Vertical  -  7W/U. 
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j)  U-component  -  8U. 

Figure  34.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  50  mph  wind  (Newport,  August  1971) (continued). 
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VflftI-9S.98  7.19 


k)    Resultant   horizontal    -    8R. 


MEAN-  «.  20     1.88 
VflRI-10.18  2.03 


mvc  m 


1)  Vertical  -  8W/U. 
Figure  34.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  50  mph  wind  (Newport,  August  1971) (continued) . 
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m)    U-component   -    9U. 


m»       nn 
HERN-55.3«2«.7V 
VRRI- 34.115  9.99 


n)  Resultant  horizontal  -  9R. 

Figure  34.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  50  mph  wind  (Newport,  August  1971) (continued) . 
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MEAN-U.U5    1.99 
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o)  Vertical  -  9W/U. 
Figure  34.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  50  mph  wind  (Newport,  August  1971) (continued) . 
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HAVE  NO 

b)  Resultant  horizontal  -  5R. 
Figure  35.   Individual  channel  horizontal  and  vertical 
wind  spectra  for  60  mph  wind  (Newport,  August  1971) . 
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c)    Vertical   -    5W/U. 
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d)  U-component  -  6U. 

Figure  35.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  60  mph  wind  (Newport,  August  1971) (continued) . 
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e)  Resultant  horizontal  -  6R. 


HPH  H/S 

MEAN=5.74     2.56 
VRRI=6.10     1.22 


"to* — I — J   I  W  Utter* — I — TTTTTWIo5 — I — TTTTtJlIo5 — J — J   4  FTTH 

HAVE  HO 

f)  Vertical  -  6W/U. 

Figure  35.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  60  mph  wind  (Newport,  August  1971) (continued) . 
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h)  Resultant  horizontal  -  7R. 

Figure  35.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  60  mph  wind  (Newport,  August  1971) (continued) . 
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j)  U-coraponent  -  8U. 

Figure  35.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  60  mph  wind  (Newport,  August  1971) (continued) . 
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k)  Resultant  horizontal  -  8R. 
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1)  Vertical  -  8W/U. 

Figure  35.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  60  mph  wind  (Newport,  August  1971) (continued) . 

96 


*fc    y',\."^'T(JLHkiil. 


"1 


T 


-+■ 


HPH  H/S 

MERN=55.96  25.01 
VRRI-UU.26  8.84 


-to* — i    i  i  nrnw — *    *  i  *nma* — n-nrrnr? — i    i  i  rnn 

HAVE  N0 

m)  U-component  -  9U. 
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n)  Resultant  horizontal  -  9R. 

Figure  35.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  60  mph  wind  (Newport,  August  1971) (continued) . 
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o)  Vertical  -  9W/U. 
Figure  35.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  60  mph  wind  (Newport,  August  1971) (continued) . 
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APPENDIX  D 
INDIVIDUAL  CHANNEL  WIND  SPECTRA:   NEWPORT,  JANUARY  1972 
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b)   Resultant  horizontal  -  2R. 

Figure  36.   Individual  channel  horizontal  and  vertical 
wind  spectra  for  40  mph  wind  (Newport,  January  1972) . 
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c)  Vertical  -  2W/U. 
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d)  U-component  -  3U. 

Figure  36.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  40  mph  wind  (Newport,  January  1972) (continued) . 
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e)  Resultant  horizontal  -  3R. 


MPH  M/S 

MEflN=ii.l8     1.87 
VARI  =  U.61      0.92 


1 — 4    A  4  iUHtr* 5 — mTTUW      T~ TTTTTUW'      I     TTTTTTr 

WflVE    NO 

f)  Vertical  -  3W/U. 
Figure  36.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  40  mph  wind  (Newport,  January  1972) (continued) . 
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h)  Resultant  horizontal  -  5R. 

Figure  36.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  40  mph  wind  (Newport,  January  1972) (continued) . 
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j)  U-component  -  6U. 

Figure  36.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  40  mph  wind  (Newport,  January  1972)  (continued)  . 
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k)    Resultant  horizontal   -    6R. 
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1)  Vertical  -  6W/U. 

Figure  36.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  40  mph  wind  (Newport,  January  1972) (continued) . 
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m)  U-component  -  7U. 
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n)     Resultant   horizontal    -    7R. 

Figure    36.       Individual    channel   horizontal   and   vertical   wind 
spectra    for    40  mph  wind    (Newport,    January    1972) (continued) . 
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o)  Vertical  -  7W/U. 

Figure  36.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  40  mph  wind  (Newport,  January  1972) (continued) . 
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b)    Resultant   horizontal    -    2R. 

Figure    37.       Individual    channel   horizontal   and  vertical 
wind    spectra   for    50   mph  wind    (Newport,    January   197  2) . 
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d)  U-component  -  3U.  \  , 

Figure  37.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  50  mph  wind  (Newport,  January  1972) (continued) . 
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f)  Vertical  -  3W/U. 

Figure  37.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  50  mph  wind  (Newport,  January  1972) (continued) . 
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h)  Resultant  horizontal  -  5R. 

Figure  37.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  50  mph  wind  (Newport,  January  1972) (continued) . 
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j)  U-component  -  6U. 

Figure  37.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  50  mph  wind  (Newport,  January  197  2) (continued) . 
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k)  Resultant  horizontal  -  6R. 
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1)  Vertical  -  6W/U. 

Figure  37.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  50  mph  wind  (Newport,  January  1972) (continued) . 
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n)    Resultant   horizontal    -    7R. 
Figure    37.       Individual   channel   horizontal   and  vertical  wind 
spectra   for   50  mph  wind    (Newport,    January   1972) (continued) . 
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o)  Vertical  -  7W/U. 

Figure  37.   Individual  channel  horizontal  and  vertical  wind 
spectra  for  50  mph  wind  (Newport,  January  1972) (continued) . 
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APPENDIX  E 
INDIVIDUAL  CHANNEL  WIND  SPECTRA:   SITKA,  NOVEMBER  197  5 
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Figure  38.   Individual  channel  horizontal  resultant  and 
vertical  wind  spectra  for  30  mph  wind  (Sitka',  November  1975) 
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c)  Horizontal  -  3R. 
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d)  Vertical  -  3W/R . 
Figure  38.   Individual  channel  horizontal  resultant  and  vertical 
wind  spectra  for  30  mph  wind  (Sitka,  November  197  5) (continued) . 
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f)   Vertical  -  4W/R. 
Figure  38.   Individual  channel  horizontal  resultant  and  vertical 
wind  spectra  for  30  mph  wind  (Sitka,  November  1975) (continued) . 
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h)   Vertical  -  5W/R. 
Figure  38.   Individual  channel  horizontal  resultant  and  vertical 
wind  spectra  for  30  mph  wind  (Sitka,  November  1975)  (continued)  . 
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APPENDIX  F 
INDIVIDUAL  CHANNEL  WIND  SPECTRA:   SITKA,  MAY  1975 
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b)   Vertical  -  1W/R. 
Figure  39.   Individual  channel  horizontal  resultant  and 
vertical  wind  spectra  for  25  mph  wind  (Sitka,  May  1976). 
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d)  Vertical  -  2W/R. 

Figure  39.   Individual  channel  horizontal  resultant  and  vertical 
wind  spectra  for  25  mph  wind  (Sitka,  May  19760  (continued) . 
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f)   Vertical  -  3W/R. 
Figure  39.   Individual  channel  horizontal  resultant  and  vertical 
wind  spectra  for  25  mph  wind  (Sitka,  May  1976)  (continued) . 
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h)   Vertical  -  5W/R. 
Figure  39.   Individual  channel  horizontal  resultant  and  vertical 
wind  spectra  for  25  mph  wind  (Sitka,  May  1976)  (continued). 
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APPENDIX  G 
INDIVIDUAL  CHANNEL  ACCELERATION  SPECTRA:   SITKA,  DECEMBER  1975 
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b)      A2. 

Figure   40.    Individual   channel   acceleration    spectra   during 

first  bridge  motion   event. 
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Figure  40.   Individual  channel  acceleration  spectra  during 
first  bridge  motion  event  (continued) . 


128 


V. 


O.OM- 


0.O2& 


*4M?" 


£1 


^jL 


y^Dfi»  itfpW  *JV»l*fc.fc4«S 


** 


25        Em        um       ZS       Cm        Cm 


niw 


-ar^ 


A  _ .  L 


e)       A5. 


O.OM- 


0     i  i.»»M-^^ 


0.02ft. 


^ 

S* 


Slw       cm" 


Iffff  >  <|^f|H^|ilfti»i  M*mi*m+++*tkm*m+m* 


Em       £i 


Cm        35" 


A.-  M .  t 


FltM 


f)       A6. 


Figure  40.   Individual  channel  acceleration  spectra  during 
first  bridge  motion  event  (continued) . 
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h)       A8. 

Figure    40.       Individual    channel   acceleration   spectra   during 
first  bridge  motion   event    (continued) . 
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Figure  40.   Individual  channel  acceleration  spectra  during 
first  bridge  motion  event  (continued) . 
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Figure  40.   Individual  channel  acceleration  spectra  during 
first  bridge  motion  event  (continued) . 
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b)   A2. 

Figure  41 


Individual  channel  acceleration  spectra  during 
second  bridge  motion  event, 
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Figure  41.   Individual  channel  acceleration  spectra 

during  second  bridge  motion  event  (continued) . 
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f)      A6. 

Figure   41.      Individual   channel  acceleration  spectra  during 
second  bridge  motion  event    (continued) . 
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Figure  41.   Individual  channel  acceleration  spectra  during 
second  bridge  motion  event  (continued) . 
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Figure  41.   Individual  channel  acceleration  spectra  during 
second  bridge  motion  event  (continued) . 
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Figure  41.   Individual  channel  acceleration  spectra  during 
second  bridge  motion  event  (continued) . 
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b)      A2. 
Figure    42. 


Individual   channel   acceleration   spectra  during 
third  bridge  motion   event. 
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Figure  42.   Individual  channel  acceleration  spectra  during 
third  bridge  motion  event  (continued) . 
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f)       A6. 

Figure  42.   Individual  channel  acceleration  spectra  during 
third  bridge  motion  event  (continued) . 
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Figure    42.      Individual   channel   acceleration   spectra  during 
third  bridge  motion   event    (continued) . 


142 


t.l 


<pwfrrtnf*  W>    t  II  I 


H 1- 


3 

r 

5 

«s — nr 


nl 


25      IB      25      «■      «S 


i)      A9, 


0.9ML 


"  tt»    MwhftftMU   I 


HtyHF  WHMMi       »    U  ■ 


H 1 1 1 1 1 1 1 


iP 


j)       A10. 


Figure    42.       Individual   channel   acceleration    spectra  during 
third  bridge  motion   event    (continued) . 
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Figure  42.   Individual  channel  acceleration  spectra  during 
third  bridge  motion  event  (continued) . 
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